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Summary. The ADM Hamiltonian formulation of general relativity with prescribed lapse 
and shift is a weakly hyperbolic system of partial differential equations. In general weakly 
hyperbolic systems are not mathematically well posed. For well posedness, the theory should 
be reformulated so that the complete system, evolution equations plus gauge conditions, is 
(at least) strongly hyperbolic. Traditionally, reformulation has been carried out at the level of 
equations of motion. This typically destroys the variational and Hamiltonian structures of the 
theory. Here I show that one can extend the ADM formalism to (i) incorporate the gauge condi- 
tions as dynamical equations and (ii) affect the hyperbolicity of the complete system, all while 
maintaining a Hamiltonian description. The extended ADM formulation is used to obtain a 
strongly hyperbolic Hamiltonian description of Einstein's theory that is generally covariant 
under spatial diffeomorphisms and time reparametrizations, and has physical characteristics. 
The extended Hamiltonian formulation with 1+log slicing and gamma-driver shift conditions 
is weakly hyperbolic. 



1 Introduction 

This paper is dedicated to Claudio Bunster in celebration of his sixtieth birthday. In a remark- 
able body of work Claudio showed that we can view the Hamiltonian formulation of general 
relativity as fundamental. (See in particular references 1261 1141 [27l .) He considered the re- 
quirement that the sequence of spatial three-geometries evolved by the Hamiltonian should 
be interpretable as a four-dimensional spacetime. From this assertion and a few modest as- 
sumptions he was able to derive the ADM Hamiltonian 1 8 1 1 of general relativity. A number 
of deep insights into the nature of gravity and matter came from his analysis, including the 
role of gauge symmetries in electrodynamics and Yang-Mills theories and the necessity for 
all matter fields to couple to gravity. 

In general relativity we are faced with the practical problem of predicting the future evolu- 
tion of strongly gravitating systems, including interacting black holes and neutron stars. Such 
problems fall into the realm of numerical relativity. Naturally, the first attempts at numerical 
modeling in general relativity were based on the ADM Hamiltonian equations. By the early 
1990's it became clear that the ADM equations were not appropriate for numerical computa- 
tion because they are not well posed in a mathematical sense. What followed was more than a 
decade of activity in which the ADM equations were rewritten in a variety of ways. The goal 
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was to produce a well posed system of partial differential equations (PDE's) for Einstein's 
theory. One strategy for modifying the ADM equations was to add multiples of the constraints 
to the right-hand sides. Another strategy was to introduce new independent variables defined 
as combinations of metric tensor components and their spatial derivatives. This later strategy 
introduces new constraints into the system, namely, the constraints that the definitions of the 
new variables should hold for all time. 

In a practical sense, the effort to re-express the ADM equations has been successful. 
Currently there are a number of formulations of the Einstein evolution equations that appear 
to be "good enough", the most popular for numerical work being the BSSN svstem. l24l l2l 
BSSN relies on a conformal splitting of the metric and extrinsic curvature. It introduces new 
independent variables, the "conformal connection functions", defined as the trace (in its lower 
indices) of the Christoffel symbols built from the conformal metric. 

BSSN and the other "modern" formulations of Einstein's theory are very clever. But at a 
basic level, these formulations are obtained through a manipulation of the equations of motion 
as a system of PDE's. What is invariably lost is the beautiful Hamiltonian structure found 
in the ADM formulation. In this paper I present a systematic procedure that can be used to 
modify the ADM equations in an effort to obtain a good system of PDE's without losing the 
Hamiltonian structure. 

A good system of PDE's is one that is mathematically well posed. As a general rule, a 
system formulated in space without boundaries must be strongly hyperbolic to be well posed. 
If boundaries are present, an even stronger notion of hyperbolicity, symmetric hyperbolicity, 
is needed to prove well posedness. We are interested in extensions of the ADM equations that, 
like the ADM equations themselves, have first-order time and second-order space derivatives. 
It turns out that a simple prescription can be given to test for strong hyperbolicity in such 
systems of PDE's. The justification for this prescription requires a rather deep mathematical 
analysis, but the prescription itself is fairly easy to apply. In section II, I discuss hyperbolicity 
and justify the prescription for well posedness with heuristic arguments. 

Another issue that has become apparent from recent numerical work is the benefit, in prac- 
tice, of incorporating the slicing and coordinate conditions (gauge conditions) as dynamical 
equations. That is, the lapse function and shift vector are not fixed a priori but are determined 
along with the other fields through evolution equations of their own. The hyperbolicity of the 
entire system of PDE's, including the equations for the lapse and shift, must be considered. 
The issues of gauge conditions and hyperbolicity cannot be separated. 

In this paper I show that the Hamiltonian formulation of general relativity can be extended 
to (i) incorporate dynamical gauge conditions and (ii) alter the level of hyperbolicity. In sec- 
tion III the ADM formulation is enlarged by the introduction of momentum variables K and p a 
conjugate to the lapse function a and shift vector /3 a . In this way the lapse and shift become 
dynamical. The new momenta are primary constraints and they appear in the action with un- 
determined multipliers A and £2". The usual Hamiltonian and momentum constraints, ^ and 
are secondary constraints. This Hamiltonian formulation of Einstein's theory is not new 
171 1121 , and is not substantially different from the original ADM formulation — like the ADM 
formulation, it is only weakly hyperbolic. This is shown in section IV. 

The Hamiltonian formulation with dynamical lapse and shift is extended in section V by 
allowing the multipliers A and fl a to depend on the canonical variables. This has two effects. 
First, it changes the evolution equations for the lapse and shift, yielding gauge conditions that 
depend on the dynamical variables. Second, it changes the principal parts of the evolution 
equations and potentially changes the level of hyperbolicity of the system. The hyperbolic- 
ity of the extended Hamiltonian formulation is considered in section VI for a fairly general 
choice of multipliers that preserves spatial diffeomorphism invariance and time reparametriza- 
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tion invariance. When the multipliers are chosen so that the evolution equations are strongly 
hyperbolic with physical characteristics, the system is equivalent in its principal parts to the 
generalized harmonic formulation of gravity. 1 9 20, 16) It is also shown that the extended 
Hamiltonian formulation with 1+log slicing and the gamma-driver shift condition is only 
weakly hyperbolic. A few final remarks are presented in section VII. 



2 Strong Hyperbolicity for Quasilinear Hamiltonian systems 

Let q„, Pji denote pairs of canonically conjugate fields. We will consider Hamiltonian sys- 
tems for which Hamilton's equations are a quasilinear system of partial differential equations 
(PDE's). Thus we assume that the Hamiltonian H is a linear combination of terms that are at 
most quadratic in the momenta and spatial derivatives of the coordinates. More precisely, H 
should be a linear combination of termsp^py, {d a q tl ){di } q v ), Pfi{d a q v ), p^, (d a q^), and 1 
with coefficients that depend on the g's[j (Here, d a denotes the derivative with respect to the 
spatial coordinates.) One would like to show that Hamilton's equations are well posed as a 
system of PDE's. The subject of well posedness is a large, active area of research in math- 
ematics and physics. In this section I present a very pedestrian account of the subject in the 
context of Hamiltonian field theory. Much more rigorous and complete discussions can be 
found elsewhere. (See, for example, references fTTI |22l l25l l6l[23l[T9l[T0l[T5l.) 

A well posed system is one whose solutions depend continuously on the initial data. For 
a well posed system, two sets of initial data that are close to one another will evolve into 
solutions that remain close for some finite time. A system is not well posed if it supports 
modes whose growth rates increase without bound with increasing wave number. A concrete 
example is given below. 

In analyzing well posedness we are primary concerned with the evolution in time of high 
wave number (short wavelength) perturbations of the initial data. For this purpose we can ap- 
proximate the quasilinear system of PDE's by linearizing about a solution. That is, we expand 
the Hamiltonian to quadratic order in perturbations, which we denote Sq^ , Sp^. We then look 
for Fourier modes of the form Sq^ = qf 1 e"°' +llc ' ,x ° /(i\k\), 8p^ = p„e"° t+ika with nonzero 
wave number k a . Here, \k\ = \J h ab k a ki, is the norm of k a defined in terms of a convenient 
metric h ab (which could be the inverse of the physical spatial metric). If the ansatz for Sq^, 
8pn is substituted into the linearized Hamilton's equations, the system becomes 

00v=(\k\A-iB-C/\k\)v (1) 

where v is the column vector v = {q\,q~2i-- -iPuPli ••-) of Fourier coefficients. Equation (0 
shows that the problem of finding perturbative modes with wave number k a is equivalent to the 
eigenvalue problem for the matrix (\k\A — IB — C/\k\). The eigenvector is v and the eigenvalue 
is the frequency a>. 

What one is really doing in the construction above is replacing the system of PDE's with 
a pseudo-differential system. The factor of i\k\ in the denominator of Sq^ is, in effect, equiv- 
alent to a change of variables in which q^ is replaced by q^j \J h ab d a d[,. In this way the 
second order (in space derivatives) PDE's are replaced with an equivalent first order pseudo- 
differential system. 1 191 

1 Note that terms proportional to d a pp_ are also allowed in H since they are related to terms 
of the form Pfi(d a q v ) through integration by parts. 



4 J. David Brown 

The behavior of CO as \k\ becomes large depends on the leading order term A in the matrix 
(\k\A — iB — C/\k\). The term A is the "principal symbol" of the system. It is constructed 
from the coefficients of the highest "weight" terms in the Hamiltonian, namely, the terms 
proportional to PfiPv, {d a qp){dbqv), Pii(d a qv) and (d a pfi). Note that it is not necessary to 
linearize the equations of motion (or expand the Hamiltonian to quadratic order) in order to 
find A. In practice we don't actually linearize, we simply identify the coefficients of the highest 
weight terms in the PDE's to form the matrix A. 

If A has real eigenvalues and a complete set of eigenvectors that have smooth dependence 
on the unit vector n a = k a /\k\, the system is said to be strongly hyperbolic. If A has real 
eigenvalues but the eigenvectors are not complete, the system is said to be weakly hyperbolic. 
It can be proved that a strongly hyperbolic system of quasilinear PDE's is well-posed. 1 191 

Here is the rough idea. Let 5 denote the matrix whose rows are the left eigenvectors of A. 
Assuming strong hyperbolicity, the eigenvectors are complete and S -1 exists. The eigenvalue 
problem {j} can be written as cov= (\k\A — iB — C/\k\)v where v = Sv, A = SAS~ 1 , B = SBS~ 1 , 
and C = SCS~ l . Note that A is diagonal with entries equal to the (real) eigenvalues. Let a 
dagger (f) denote the Hermitian conjugate (complex conjugate * plus transpose T). Since 
A 1 " = A, we find 

(CO — £0*)v T v = v 1 (o)v) — (cov)^v 

= v ] '(\k\A-iB-C/\k\)v-v ] '(\k\A + iB T -C T /\k\)v 

= -itf{M+M*)v (2) 

where M = B — iC/\k\ . The left-hand side includes the factor (co — CO*) = 2iSco = — 2i^i(ico), 
so Eq. (|2j can be written as 29t(/co) = v^M+M^v/v^v. It follows that 9t(z'(u) < where 
2t~ 1 is the maximum over \k\ of the matrix norm of (M + M^). [The matrix norm is the 
maximum of the real number v^(M+M^)v/(v^v).] From this argument we see that the growth 
rate for the mode k a is bounded; it can grow no faster than e'/ T where T is independent of k a . 

Consider a simple example in one spatial dimension with two pairs of canonically conju- 
gate fields, q\ , p\ and q2, P2- Let the Hamiltonian be given by 

H = J dx^[( Pl f + (P2) 2 + (q' l ) 2 + (q' 2 ) 2 }+2p 2 q' l +2p 1 q' 2 

+Pl(<72 + 9l)+P2(92-<?l)+<?l<72 j ■ 

Hamilton's equations are 

q\ = Pi +2</ 2 +<?2+<7l , 
42 = P2 + 2<7i +qi-qi , 

pi = q'[ + V2 + pi - pi - 91 , 
Pi = q'{ + Vi -pi- pi - q\ ■ 

In this example the PDE's are linear so the linearization step is trivial: q^ -4 Sq^, p^ 
Now insert the ansatz 8q^ = q ft e ico '+ ikx /(i\k\), 8 Pll = p tl e" s >' +ikx . This yields 

(x)q\ = \k\(2nq 2 +Pl)-i{qi+q2) , 
(Oq 2 = \k\(2nqi + p 2 ) - i{q2 - qi) , 
®Pl = \k\(qi+2np 2 )-i(p 2 - pi)+q 2 /\k\ , 
(°P2 = \k\{q2 + 2npi)+i(pi +p 2 )+q l /\k\ , 



(3) 



(4a) 
(4b) 
(4c) 
(4d) 

S Pu- 
pa) 
(5b) 
(5c) 
(5d) 
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where n = k/\k\ is the sign of the wave number k. Collecting the unknowns into a column 
vector v = (q 1 ,qi,p\ ,Pi) T , we see that these equations become cov = (\k\A — iB — C/\k\)v 
where the matrices are given by 



/ 2n 1 \ 
2n 1 
1 In 
\ 1 2n j 



1 1 \ 
-110 
0-11 
0-1-1/ 



(° 



V-i 



0\ 

-10 
0/ 



(6) 



(7) 



The principal symbol A has real eigenvalues ±1, ±3, and a complete set of eigenvectors. 
Therefore this system is strongly hyperbolic. The modes with wave number k have frequencies 
co = ±\k\ + ^(l/|fc|) and co = ±3\k\ + G(l/\k\). In particular the imaginary parts of the co's 
do not grow with increasing \k\. 

Now replace the terms 2p2<?'i + 2piq' 2 in the Hamiltonian with pi4\ ~ P\4i- The principal 
symbol becomes 

fO-n 1 0\ 
n 1 
1 n 
\0 1 -n 0/ 

while B and C are unchanged. The eigenvalues of A vanish and there are only two independent 
eigenvectors. Therefore this system is weakly hyperbolic. The modes with wave number k 
have frequencies co = ±i^/2\k\ + Gil) and co = ±^/2\k\ + The modes with frequency 
co ~ — iyj2\k\ will grow in time at a rate that is unbounded as \k\ increases. 

The system described by this last example is not well posed. Indeed, consider two ini- 
tial data sets that differ from one another by terms ~ q^e ,kx /\k\ 2 , ~ pfie' kx /\k\ where 
(?l > 92 j Pi, Pl) T is an eigenvector with eigenvalue co = —iy/2\k\ + 0(1). In the limit as 
\k\ — ► 00 these terms vanish and the two initial data sets coincide. However, if we evolve 

these data sets the solutions will differ at finite time t by terms q^ ~ q li e^^ t+lkx /\k\ 2 , 

Pii ~ pjie^^ t+lkx /\k\. These terms do not vanish in the limit \k\ — » 00. This system is ill 
posed because the solution at finite time does not depend continuously on the initial data. 

In some cases it may be possible to model a physical system with ill posed PDE's and 
to gain important physical insights through a formal analysis. Claudio's beautiful work on the 
(weakly hyperbolic) ADM equations is a perfect example! One can imagine that the initial data 
are analytic, in which case the Cauchy-Kowalewski theorem guarantees that a solution exists 
for a finite time. But most data, even smooth data, are not analytic. From a computational point 
of view, having an ill posed system is unacceptable. Numerical errors will always introduce 
modes with large wave numbers, with the size of \k\ limited only by the details of the numerical 
implementation. For example, with a finite difference algorithm the maximum \k\ is roughly 
1 j Ax where Ax is the grid spacing. In practice it does not take long for the numerical solution 
to become dominated by this highest-wave number mode. As the grid resolution is increased 
(Ax is decreased), the unwanted highest wave number mode grows even more quickly. For 
practical numerical studies, we need our system of PDE's to be well posed. 

The analysis outlined above leads to the following test for strong hyperbolicity. We begin 
by constructing the principal symbol A from the principal parts of Hamilton's equations. The 
principal parts of the q^ equations are the terms proportional to p^ and d a q^ . In these terms we 
make the replacements p^ — ► p^ and d a q^ — » n a q^. The principal parts of the p^ equations 
are the terms proportional to d a p^ and dadblp.- in these terms we make the replacements 
d a pfx — > n a pn and d a db1ii ~ * n anbqfi - The principal symbol A is the matrix formed from the 
coefficients of the if s and /5's. The next step is to compute the eigenvalues and eigenvectors 
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of A. If A has real eigenvalues and a complete set of eigenvectors that depend smoothly on n a , 
the system is strongly hyperbolic and the initial value problem is well posed. 



3 ADM with dynamical lapse and shift 

The Einstein-Hilbert action is S = / d x^/^gR where g is the determinant of the spacetime 
metric and R is the spacetime curvature scalar. Units are chosen such that Newton's constant 
equals 1 / ( 16tt) . With the familiar splitting of the spacetime metric into the spatial metric g ab , 
lapse function a, and shift vector /3 a , the action becomes 

S[g,a,P} = Jd 4 xa^g-(R + K ab K ab -K 2 ^ . (8) 

The extrinsic curvature is defined by 



l -[g ab -2D {a l5 b) ) , (9) 



Kab = -2a 

and D a denotes the spatial covariant derivative. The Hamiltonian can be derived in a straight- 
forward fashion if one recognizes that the action does not contain time derivatives of the lapse 
and shift. Time derivatives of the spatial metric appear through the combination K ab . Thus, we 
introduce the momentum 

P ab = ^-=Vg(Kg ab -K ab ) , (10) 



ah 

where the Lagrangian density Jzf is the integrand of the action. This definition can be inverted 
for g ab as a function of P ab and used to define the Hamiltonian: H = J d 3 x (P ab g ab — Jif) . 
This yields the ADM Hamiltonian 

H = J d i x{aj^' + fi a ^ a ) , (11) 

where 

= -L {f* h P ab - P 2 /l) - JgR , (12a) 

M a = -2D b P b a (12b) 

are the Hamiltonian and momentum constraints. 

In the analysis above the lapse and shift are treated as non dynamical fields. They appear 
in the Hamiltonian form of the action, 



%, P, a, j8] = f t ' f dt j d\ [pobgah - a.J^ - j8 a ^ } 



(13) 



as undetermined multipliers. Here f,- and tt are the initial and final times. Extremization of S 
with respect to a and j5 a yields the constraints ,yf = and ./£ a = 0. 

We can enlarge the ADM formulation to include the lapse and shift as dynamical 
variables. f7lll2l Consider again the action In addition to the momentum P ah conjugate 
to the spatial metric, we also define conjugate variables for the lapse and shift: 
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7T=^=0, (14a) 
da 

d Sf 

Pa = — ^=0. (14b) 

This leads to primary constraints n = and p fl = 0. The resulting Hamiltonian is not unique; 
it is only determined to within the addition of arbitrary multiples of the constraints: 

H = Jd 3 x{aJff + fi a Jta+ATZ + £l a p a ) . (15) 

The coefficients A and Q a appear as undetermined multipliers in the action, which now reads 
S\g,P,a,x,p,p,A,Q] 

= jf' J ' dt J d 3 x^P ah g ab + xa + p a p a -ajr-p a ^ a -AK-n a p a } . (16) 

The equations of motion, 55 = 0, ar^3 

Sab = Sfpgab + ^= (2Pab ~ Pgab) , (17a) 

P ab = ^ p P ab + (<5 t fl S* - f b g cd /A) (PP cflf - 2/> ef P<' ) 

-a^gG ab + ^(D a D h a-g ab D c D c a) , (17b) 

& = A , (17c) 

* = -Jt? , (17d) 

j3" = i2 fl , (17e) 

p a = -J( a , (170 

» = , (17g) 

P„ = , (17h) 

where G ai denotes the spatial Einstein tensor and Jfa is the Lie derivative with respect to J3°. 

The equations above must hold for each time t; <t <tf. They are equivalent to the Ein- 
stein equations supplemented with evolution equations for the lapse function and shift vector. 
In particular, observe that 71 and p a must vanish for all time by Eqs. ( |17g|17ht . It follows 
that the time derivatives of n and p a must vanish. In turn, Eqs. J17dll7ft imply that the usual 
Hamiltonian and momentum constraints are zero. Eqs. l |17a|17b| ( are the familiar ADM evo- 
lution equations, and Eqs. ( 1 17cl 1 7et supply evolution equations for the lapse and shift. 

The evolution equations <17h -f) are Hamilton's equations derived from the Hamiltonian 
d!5t . The time derivative of any function F of the canonical variables is F = {F.H} where the 

fundamental Poisson brackets relations are defined by {gab{ x )^P cd { x ')} = ^(a^j) 3 (•*•>•*')> 
{a{x),K(d)} = S 3 (x,x > ), and {f5 a (x),p b (x')} = 8^8 3 (x,x l ). We can interpret Hamilton's 
equations as an initial value problem by following Dirac's reasoning for constrained Hamil- 
tonian systems. 1 13] The initial data are chosen such that the primary constraints 71 and p a 
vanish at the initial time. According to Eqs. J17dll7ft . these constraints will remain zero as 
long as 34? and ^ a are constrained to vanish as well. Thus we impose 3if = and = as 
secondary constraints. The complete set of constraints, % = 0, p a = 0, = 0, and d( a = is 
first class. 



Throughout this paper I ignore the issues that arise when space has boundaries.! 21 3 
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4 Hyperbolicity of ADM with dynamical lapse and shift 

Hamilton's equations dl7h-f) are equivalent to the ADM equations plus evolution equations 
a — A, jS a = £i a for the lapse and shift. Let us analyze the level of hyperbolicity of these 
PDE's. The principal parts of the q equations are the terms that are proportional to p's or 
first spatial derivatives of q's. The principal parts of the p equations are the terms that are 
proportional to first spatial derivatives of p's or second spatial derivatives of q's. Thus, we find 



k)gab = 2g c{a d h )P C + -^(ZPab-Pgab) , (18a) 

+^-s ab s cd g ef {dcd egdf - d c d dgef ) 

+Vg(g ac g M -8 ah g cd )dcd d a, (18b) 

d a ^ -P a d d a , (18c) 

do* = Vgg ah g Cd (dad c gbd - d a d hgcd ) - p a d a X , (18d) 

d^ a s -p b d b $ a , (18e) 

d pa = 2gacd b P hC - P b d b p a , (18f) 



where the symbol = is used to denote equality up to lower order (non principal) terms. These 
equations have been expressed in terms of the operator do = d r — fi"d a so that the characteristic 
speeds (the eigenvalues of the principal symbol) are defined with respect to observers who are 
at rest in the spacelike slices. 

We now construct the eigenvalue problem fiv = Av for the principal symbol A. The prin- 
cipal symbol is found from the coefficients on the right-hand sides of Eqs. d!8t . These coeffi- 
cients are divided by a factor of a so that the characteristic speeds will be expressed in terms 



of proper time rather than coordinate time. The result is 

Vgab = -n (a fo) + ^p[2Pab-gab(Pnn+PABS AB )] , (19a) 

UPab = ^[gab - 2n {a g b)n +n a n b {g nn + g AB S AB ) - g ab gAB& AB ] 

-^{gab-n a n b )a , (19b) 

= -(/} -n/a)& , (19c) 

lilt = -^-g AB S AB -(p -n/a)Tt, (19d) 

HPa = -(P-n/a)p a , (19e) 

Wa = -Pna-(P-n/a)p a . (19f) 



where fl is the eigenvalue and v = {g ab ,P ab ,cc, KiPa,Pa) I s me eigenvector. In these equations 
n a is normalized with respect to the spatial metric, n a g ab n b = 1, and a subscript n denotes 
contraction with n a . We have also introduced an orthonormal diad e a A with A = 1,2 in the 
subspace orthogonal to n a . That is, n a e A = and e^g ab e b B = 5 A b- A subscript A on a tensor 
(such as the metric g ab or momentum P ab ) denotes contraction with e a A . 
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The eigenvalue problem d 1 9b splits into scalar, vector and trace-free tensor blocks with 
respect to rotations about the normal direction n a . The scalar block is 

Mln« = -A. + 4= ( P »» - P AB$ AB ) , (20a) 
a \/8 

llgABS™ = -^pPnn , (20b) 

llPnn = , (20c) 

LiPabS ab = ~ Vssab8 ab ~^Vsa, (20d) 
jia = — (/3 • n/a)a , (20e) 

Lin = -^-g AB 8 AB -(p-n/a)K, (20f) 

= -<fi-n/a)fr, , (20g) 
2 - 

MP« = - (j8 • «/a)p„ . (20h) 

This block has eigenvalues and — (]3 -n/a), each with multiplicity 4. There is only one 
eigenvector with eigenvalue 0, so the eigenvectors are not complete. The vector block is 

1 - 2 

MfnA = -Pa + —pP„A , (21a) 

M^A = , (21b) 

llp~ A = -(^.n/o)fii, (21c) 

2 . 

MPa = -^«a - (0 • »/«)Pa ■ (21d) 

It has eigenvalues and — (J3 • each with multiplicity 2. There is only one eigenvector 

with eigenvalue 0, so the eigenvectors are not complete. Finally, the trace-free tensor block is 

PIab = ^fL - (22a) 

A = ■ (22b) 

This block has eigenvalues ±1 and a complete set of eigenvectors. Because the eigenvectors 
for the scalar and vector blocks are not complete, the system (T7J is weakly hyperbolic. 



5 Extending the ADM formulation 

In the previous section we modified the ADM Hamiltonian formulation of general relativity 
so that the lapse function a and shift vector j8 fl are treated as dynamical variables. Their 
canonical conjugates are denoted n and p a , respectively. The undetermined multipliers for the 
constraints n = 0, p a = are A and Q a . These multipliers are freely specifiable functions 
of space and time. They determine the slicing and spatial coordinate conditions through the 
equations of motion a — A and p" = £1". 

Here is the key observation. The histories that extremize the action are unchanged if we 
replace the multipliers by A — * A + A and Q a — > Q a + £2 a , where A and £l a are quasilin- 
ear functions of the canonical variables. By quasilinear, I mean that the principal parts of A 
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and Q" are linear in the momenta (P ab , K and p a ) and first spatial derivatives of the coor- 
dinates (d c g a j,, d c a and (3 c j8 fl ) with coefficients that depend on the coordinates. With these 
replacements the action becomes 



S\g,P,a,n,p,p,A,£i] = J dt J d i x^P ab g ab + na + p a li a -aJf-l5 a .^' a 

-(A+A)x-(fl a + £2 a )p a ^ , (23) 

and the Hamiltonian is 

H= [d 3 x(a,J^ + p a ^ a + {A+A)7t + (Q a + Q a )p a ) . (24) 



The solutions to the equations of motion are unaltered because the extra terms are proportional 
to the constraints % — 0, p a = 0. In the Hamiltonian formulation we can dispense with the 
original multipliers A and Q. a \ that is, these quantities can be absorbed into the functions A 
andiK 

With A and Q a restricted to be quasilinear in the momenta and first spatial derivatives of 
the coordinates, the equations of motion become 

a . dA dQ c 

gab = Jffigab + -JZ (2Pab - Pgah) + * + J^Pc , (25a) 

pab = j? pP ab + a {§ a 5 b_ g ab gcd/4){pp cd_ 2p ce p d ) 
\Jg 

-a^g~G ab + s /g'(D a D b a - g ab D c D c a) 

dA dfr _ / dA \ ^ ( dfr \ /ncux 
-j. — n- 3 — p c + d d — -TZ \+d d \ — -p c , (25b) 

dgab dgab \d(d d g ab ) J \d(d d g ah ) J 

a = A+A + —K+^—p c , (25c) 
.„ dA d& n / dA \ ^ ( dfr \ 

* = -^-^-^P^ 9 " 1 2(1 „\ n ) +d d\ J?3 ^ Pc . ( 25d ) 

da da \d(d d a) I \d(d d a) I 

P a = Q a + n a + ^-n + ^-p c , (25e) 

dp a dp a 

„ dA dfr . / dA \ ^ ( dfc \ 

^ = -^-df^~W Pc + dd {WaW) n ) +dd {WaW) Pc ) ' (25f> 

K = , (25g) 

Pa = , (25h) 

Equations J25al25bb are the usual ADM equations apart from terms proportional to the con- 
straints, n = and p a = 0. The equations that govern the slicing and spatial coordinates are 
generalized by the presence of the functions A and Q a . Apart from terms that vanish with 
the constraints n = 0, p a = 0, the evolution equation for the lapse becomes a = A + A and 
the evolution equation for the shift becomes p a = Q" + £2 a . The equations for it and p a are 
modified, but once again we see that the complete set of constraints % = 0, p a — 0, M 1 = 0, 
and = is first class. 
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In principle we can choose A and £2 a to be any set of quasilinear functions of the canonical 
variables. In practice we might want A and £2" to satisfy certain transformation properties. 
For example we can restrict A to be a scalar and £2 a to be a contravariant vector under spatial 
diffeomorphisms. This allows us to maintain a geometrical interpretation of the equations of 
motion. In particular this allows us to prepare and evolve identical geometrical data using 
different spatial coordinate systems. 

Another property that can be imposed on the formalism is reparametrization invariance. |13l 
This is invariance under a change of coordinate labels / for the constant time slices. Consider 
the infinitesimal transformation t — > t — e(t), In the usual ADM system, the variables g a i, and 
P ab transform as scalars under time reparamerization: 8g a [, = Eg a b and 8P° b = eP ab . The time 
derivative of the metric, g a j,, transforms as a covariant vector. In one dimension a covariant 
vector transforms in the same way as a scalar density of weight + 1 . It follows that the term 
P" b g c ,b that appears in the action is a weight + 1 scalar density. For reparametrization invari- 
ance to hold, the integrand of the action should transform as a weight + 1 scalar density since 
it is integrated over t. In particular the lapse function a and shift vector j8 fl , which multiply 
the scalars and j$ a , must transform as scalar densities of weight +1. 

Observe that the time derivatives a and f5 a are constructed from coordinate derivatives of 
scalar densities and, as a consequence, these quantities do not transform as tensors or tensor 
densities. We need to replace the coordinate derivatives (dots) with covariant derivatives. We 
can do this by choosing a background metric for the time direction. This should be viewed 
as part of the gauge fixing process. Now, the physical metric for the time manifold is a 2 , so 
let a 2 denote the background metric. The covariant derivative built from a 2 , acting on the 
densities a and is defined by 

a = a — (a/a)a, (26a) 
/3 fl = $ a -{a/a)P a . (26b) 

The extra terms needed for reparametrization invariance can be built into the action or Hamil- 
tonian by including a term (St /a) a in the function A and a term (a/a)f} a in the function 
£2". 

With the appropriate terms included in A and £2" , the time derivatives of the lapse and 
shift appear in the action only in the combinations nh and p a $ a - Since a and jS° are covariant 
vector densities of weight +1, we see that 71 and p a must transform as contravariant vectors 
with no density weight. In one dimension, contravariant vectors transform in the same way as 
a scalar density of weight — 1 . We will consider n and p a to be scalar densities of weight — 1 
under time reparametrization. It follows that, apart from the terms (a /a) a and (a/a)p a , the 
multipliers A+A and £2" + £2" should transform as scalar densities of weight +2. 

We have now established the rules for adding terms to the functions A and £2" such that the 
resulting formulation is invariant under time reparametrizations: these terms must be weight 
+2 densities built from the scalars g a \,, P ab , the weight +1 densities a, j8 a , and the weight 
— 1 densities %, p a . We can also insist that A and £2° should be, respectively, a scalar and 
a contravariant vector under spatial diffeomorphisms. With these properties in mind, a fairly 
general form for A is 

2 3 

A = (&./a)a + fi a D a a-Ci—P + C A — iz. (27) 
s/g 

The first term in required for reparametrization invariance. The second term will allow us 
to combine shift vector terms into a Lie derivative Jzfjrj acting on a. The terms multiplied 
by constants C\ and C4 are principal terms that will affect the hyperbolicity of the resulting 
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equations. There are other principal terms that one can add, such as a 2 j8 fl p a j \fg, but the form 
above will be general enough for present purposes. There are also lower order terms that one 
can add to A . 

For £2" we must construct a spatial vector that [apart from the term (a/a)fi a ] transforms 
as a weight +2 density under time reparametrizations. There are several obvious ways to 
construct a spatial vector from the canonical variables at hand. There are some possibilities 
that are not so obvious. Recall that the difference of two connections is a tensor. Thus, the 
combination T b — fL" is a spatial tensor if E" are the Christoffel symbols built from the 
physical metric g ab and T bc are the Christoffel symbols built from a background metric g ab . 
In setting up a numerical calculation, for example, on a logically Cartesian grid, we can take 
g ab to be the flat metric in Cartesian coordinates. Again, we view the introduction of the 
background structure g ab as part of the gauge fixing process. 

The general form for tl a that we will consider is 

Q a = (&/a)p a +p h D h l3 a +C 2 a 2 (r b " c -f b a c )g bc 

+C 3 a 2 (r c l-ff b )g ab -C 5 aD a a-C 6 ^-p a . (28) 

where D a is the covariant derivative compatible with g ab . The first term is required for time 
reparametrization invariance. The second term will allow us to combine time derivatives and 
shift vector terms into the operator do in the principle parts of the equations for j3 a and p a . 
The remaining terms will modify the principal parts of the equations of motion and can affect 
the hyperbolicity of the system. There are other principal terms that we could add to Q", such 
as a 2 nP a /y/goi aP ab p b /^g. We can also add lower order terms. 

With these expressions for A and ti a , we find the following equations of motion by vary- 
ing the action ( I231 >: 

2 

(X cc 

gab = -2/3 gab + — | (2P a b - Pgab )-C l —Kg ab , (29a) 

P c,b = ^pP° b + {%S b d - g ah g cd /4) (PP cd - 2P ce P d e ) - a^G ab + 

^fg{DPD b a - g ab D c D c a) +Ci^= (P ab - Pg ab /2)x + C 4 ^- K 2 g' 

Vg 

+C 2 D {a [p b )a 2 )-C 5 ap (a D^ a + X - (C 3 - C 2 )D C (a 2 p c )g ab 
+c 2 a 2 ft K-S)rs M +C3« 2 (rj-i:3)piV )c 

-C b ^(p a p b + p c p c g ab /2), (29b) 

2 3 

a = 3f B a+A -C\— P + 2C 4 — n . (29c) 

2 

n = J2f p % - MP + 2Ci —Pn - 3C 4 — n 2 - 2C 2 a {Ef c - f b a c )g bc p a 



1 r ,ab 



-2C 3 a(r b b -f b b )p a -C 5 aD a p a + 3C 6 ^-p a p a , (29d) 

\Jg 



2 
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-C 5 aD a a-2C 6 — p a , (29e) 

Vs 

p a = D b (p b p a ) - p h D a p b - M a - nD a a , (29f) 
n = , (29g) 

Pa = . (29h) 

These equations are generally covariant under spatial diffeomorphisms and time reparametriza- 
tions. Equations J29h-f) are generated through the Poisson brackets by the Hamiltonian i24\ . 
Note that for g a b and P ab , which are scalars under time reparametrization, the covariant time 
derivative (circle) is equivalent to a coordinate time derivative (dot). 



6 Hyperbolicity of the Extended ADM formulation 

The principal parts of the extended ADM evolution equations J29b-f) are: 

2 

d gab = 2g c{a d b) p c + ^- (2P ab - Pg ah ) - C x ^—ngab , (30a) 
V 8 \f& 

^pab c* ^-g^g bd g^\d e d f g cd -2d e d {c g d)f + d c d dgef ) 

+ ^8 ab 8 cd g ef (dcd e g df - d c d dgef ) + Vg(g aC S hd - g ah 8 cd )dcd d a 

+a 2 [c 2 g c<a g h)d + (C 3 - C 2 )g ab g cd /2\ d cPd , (30b) 

2 3 

d Q a = -Ci— P + 2C A — % , (30c) 
Vs Vg 

d Q x = Vgg ah g cd (dadcgbd - d a d b g cd ) - C 5 ag ab d aPh , (30d) 
d p" - a 2 [c 2 g ac g bd + (C,-C 2 )g ab g ed /2\ d hgcd 
a 3 

-C 5 ag ab d h a-2C 6 — p a , (30e) 

Vg 

k)Pa = 2g ac d b P bc , (30f) 

From here it is easy to construct the eigenvalue problem for the principal symbol. The symbol 
decomposes into scalar, vector, and trace-free tensor blocks under rotations orthogonal to the 
normal vector n a = k a /\k\. For the scalar sector, we find 

ilgm = -p„ + -^-(P m -P A BS AB )-C l -^7i , (31a) 

« Vg Vg 

UgAB 8 AB = - ^pPrm -2C\-^-K , (31b) 

Vg Vg 

\iPrm= y(C 3 +C 2 )p„ , (31c) 

HP AB 5 AB = -^V8gAB8 AB -^a+a{C 3 -C 2 )p„ , (31d) 

2 

Ai« = -C l -^ = (P m +P AB 8 AB ) + 2C 4 ^-7t, (31e) 

Vg Vg 
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fi7t = -^-g AB 8 AB -C 5 p n , (31f) 

2 

\x.p\ = ^(C i +C2)gnn + ^(C3-C 2 )gABS AB -C 5 a-2C 6 ^-p„ , (31g) 
2 2 ^/g 

MP* = ^Pnn ■ (31h) 

Again, the subscripts n and A denote contraction with n" and e^, respectively. The vector block 
is 

1 - 2 

Hg nA = -P A + —P nA , (32a) 



a 

and the trace-free tensor block is unmodified from before 



f-P«A = T C 2PA , 02b) 



= C 2 ag„A-2C 6 —p A , (32c) 

Vs 

2 - 

PPa = -PnA , (32d) 



Hi/ B = -j=F& , (33a) 
A = ^fm ■ 03b) 



The eigenvalues for the scalar block are ±^2C\ and ±\/C 2 +C3, each with multiplicity two. 
The eigenvalues for the vector block are ±\/C5, each with multiplicity two. For the tensor 
block the eigenvalues are ±1 and the eigenvectors are complete. 

The eigenvalues are the characteristic speeds with respect to observers at rest in the space- 
like hypersurfaces. Let us choose the C's such that the characteristics are ±1, that is, along 
the physical light cone. Then we must have C\ = 1/2, C 2 = 1, and C3 = 0. With this choice a 
careful analysis of the scalar block shows that the eigenvectors are complete only if C4 = 1/8, 
C5 = 1, and C(, = 1/2. For these values of the constants the eigenvectors for the vector 
block are complete as well. It can be shown that with these values for the C's the principal 
parts of the equations of motion are equivalent to the generalized harmonic formulation of 
relativity. [9 20 1 6 1 This will be discussed elsewhere. [ 5 1 

The Hamiltonian system l |29t is strongly hyperbolic with the choice of C's above. For this 
formulation the gauge conditions are 

a 2 P a 3 K 

& = p a D a a+A-—— + ——, (34a) 

2 4 y/g 

p a = p h D b p a +Q a + a 2 (r^' c -f£)g bc - aD a a - ^-p tt (34b) 

v ° 

If we choose A and £l a to vanish, these gauge equations become 

a — p a d a a = a(a/a) — a 2 K , (35a) 
$ a -p b d h p a = p a (&/a) + a 2 (r b a c -r h a c )g bc -ag ab d b a , (35b) 
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to within terms that vanish with the constraints it = 0, p fl = 0. Here, K = P/(2y/g) is the trace 
of the extrinsic curvature. These gauge conditions agree, to within lower order (non-principal) 
terms, with the gauge conditions for the generalized harmonic formulation of gravity. 11 61 

It is not difficult to find choices for the constants (the Cs) that make the system J29b 
strongly hyperbolic. The trick is to find a strongly hyperbolic system with desirable coor- 
dinate and slicing conditions. Conditions similar to l !35t . supplemented with "gauge driver" 
equations, have been applied to the binary black hole problem with mixed results. 1 20 . .161 [171 
With the BSSN evolution system, the gauge conditions that work well for the binary black 
hole problem are "1+log slicing", 

a-p a d„a = -2aK , (36) 

and "gamma-driver shift". The gamma-driver condition is usually written as a system of two 
first order (in time) equations for /3" and an auxiliary field B a = 4j8 fl /3. These equations, along 
with suitable initial conditions, can be integrated to yield a single first order equation for the 
shift. [28 1 Expressed as either two equations or one, the gamma-driver condition depends on 
the trace (in its lower indices) of the Christoffel symbols built from the conformal metric. In 
terms of the physical metric, we can write the single-equation form of the gamma-driver shift 
condition as 

1 

3 i 

where tj is a constant parameter. The first term on the right-hand side (apart from the factor 
of 3/4) is the trace of the conformal Christoffel symbols. 

Let us see if we can find a set of functions A and £2" that yield the gauge conditions above, 
and then ask whether the resulting system is strongly hyperbolic. In this example we dispense 
with any attempt to construct a formulation that is covariant under spatial diffeomorphisms or 
time reparametrizations. Comparing Eqs. d25c|25eb with the 1+log slicing condition l |36t and 
the gamma-driver shift condition l |37| l, we find 

A = p^a-aP/^/g+Fm+Ffpa , (38a) 



P a -P%P a = ly/f /3 



n a c s hc +-,s a % c c 



nP a , (37) 



£2 a = P»d h p a + \jg 2li 



r b a c g hc + lg ab r b c c 



Tlp a +Ffx + Ffp b , (38b) 



where the F's are functions of the coordinates g ab , a, and p a (and not their spatial deriva- 
tives). In this case the analysis of hyperbolicity is complicated by the fact that the scalar and 
vector blocks of the principal symbol are coupled by the terms F^F?, and F? b . The combined 
scalar/vector block has real eigenvalues, but completeness of the eigenvectors can be achieved 
only if F% h depends on the normal direction n a = k a /\k\. This is not acceptable; the Hamilto- 
nian cannot depend on the propagation direction of a perturbative solution. We conclude that 
there is no Hamiltonian of the form ( 124b that yields strongly hyperbolic equations with 1+log 
slicing d36t and gamma-driver shift l !37t . 



7 Concluding Remarks 

This paper outlines a procedure for constructing Hamiltonian formulations of Einstein's the- 
ory with dynamical gauge conditions and varying levels of hyperbolicity. One can use this 
procedure as a tool to help identify well posed formulations of the evolution equations that 
also maintain Hamiltonian and variational structures. The issues of dynamical gauge condi- 
tions and hyperbolicity cannot be separated. They are both dictated by the dependence of the 
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multipliers A and Q a on the canonical variables. There are many possibilities that one can 
explore for this dependence. 

The Hamiltonian and variational formulations of general relativity have shaped our per- 
spective and provided deep insights into the theory. In addition, there are a number of practical 
uses for a Hamiltonian/variational formulation. With an action principle we can pass between 
spacetime and space-plus-time formulations by adding or removing momentum variables. We 
can develop a fully first order multisymplectic version of the theory. We can also develop new 
computational techniques such as variational and symplectic integrators. 11 81 |4l 

One important issue that has not been addressed here is the constraint evolution system. 
In numerical simulations it is important to control the growth of constraint violations. This 
might be accomplished in the present framework by including appropriate terms in A and Q a 
to ensure that the constraints are damped. For example, a damping term — Cn (where C is a 
positive constant) can be added to the it equation by including a lower order term Ca in A . 
This issue will be explored in more detail elsewhere. 1 5 1 

The formalism outlined here can be further extend by introducing dynamical equations for 
A and Q a . This is accomplished by introducing momentum variables conjugate to these multi- 
pliers. The new momentum variables are primary constraints and are accompanied by a new set 
of undetermined multipliers. Dynamical equations for A and £2" are introduced by allowing 
the new multipliers to depend on the canonical variables. In this way we can construct gauge 
driver conditions similar to the ones used with the generalized harmonic formulation. 1 2011X71 
We can also allow for gauge conditions that are expressed as systems of PDE's, such as the 
two-equation version of the gamma-driver shift condition. 
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